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Abstract. The alignment of the Fermi level of a metal electrode within the gap of the highest 
occupied and lowest unoccupied orbital of a molecule is a key quantity in molecular electronics. 
Depending on the type of molecule and the interface structure of the junction, it can vary the 
electron transparency of a gold/molecule/gold junction by at least one order of magnitude. In 
this article we will discuss how Fermi level alignment is related to surface structure and bonding 
configuration on the basis of density functional theory calculations for bipyridine and biphenyl 
dithiolate between gold leads. We will also relate our findings to quantum- chemical concepts 
such as electronegativity. 



1. Introduction 

Interest for electron transport in nano-scale contacts has recently intensified, because (i) the 
advent of the technologically motivated field of molecular electronics, where the use of organic 
molecules as the active parts in electronic devices is envisioned, QIISI (ii) recent progress in the 
experimental techniques for manipulating and contacting individual molecules [SlEl, and (iii) 
the availability of first principles methods to describe the electrical properties of single molecule 
junctions on an atomic scale with high accuracy. 5, 6 These latter methods are usually based 
on density functional theory (DFT) in combination with a Keldysh formalism. 

The alignment of the Fermi level of a metal electrode within the gap of the highest occupied 
(HOMO) and lowest unoccupied orbital (LUMO) of a molecule is a key quantity in molecular 
electronics, which can vary the electron transparency of a single molecule junction by an order 
of magnitude. In a recent publication 7 we presented a quantitative analysis of the relation 
between this level alignment and charge transfer for bipyridine and biphenyl dithiolate (BPDT) 
molecules attached to gold leads based on density functional theory (DFT) calculations. The aim 
of the current work is to extend this study by analyzing the effect of varying surface structure 
and bonding configuration, where we also relate our findings to the quantum-chemical concept 
of electronegativity. 

2. Method 

All electronic structure calculations in this study are performed using a plane wave 
implementation of DFT 10 with an energy cutoff of 340 eV, where the effect of the core 
on the valence electrons is approximated by ultra-soft pseudopotentials 8 , and exchange and 
correlation potentials and energies are parametrized as a PW91 functional. 9^ In our calculations 
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Figure 1. (Color online) a) Variation of the surface contact geometry for a bipyridine molecule 
attached to a Au (111) surface, where Nc (see text) varies from three to nine. Smaller dark 
spheres define the nitrogen positions, b) Transmission functions, c) energy shifts of MOl (see 
text) depending on the distance d, where do (2.18 A ) is the equlibrium bond length and d) 
charge dependent energetic positions of MOl in an isolated molecule mimicking the interaction 
with the three different structures in a). All energies are given relative to the Fermi energy of 
the Au surfaces. 

the supercells for the central region are defined by 3 x 3 atoms in the directions perpendicular 
to the transport direction, and we used a 4 x 4 k-point grid in the transverse Brillouin zone. 

The conductance of the molecular junctions is calculated using a general non-equilibrium 
Green's function formalism for phase-coherent electron transport. ^TjThe Green's function of 
the central region, containing the molecule and three to four layers of the gold electrodes, is 
evaluated in terms of a basis consisting of maximally localized Wannier functions which are 
defined by an appropriate transformation of the Kohn-Sham eigenstates. [E] The couplings to 
the semi-infinite gold electrodes are included via self-energy matrices which are also represented 
in the same Wannier function basis. 

3. Transmission functions and MO level shifts 
3.1. Bipyridine 

In a recent article [T^ we presented a theoretical analysis of the conductance of a bipyridine 
molecular junction inspired by measurements on the same system. [14 Since the atomic 
configuration of the electrode's surfaces were not known, we scanned a wide range of surface 
structures and found a strong dependence of our results on their geometries but overall agreement 
between theory and experiment in orders of magnitude (^ 0.01 Go). All our calculations have 
been performed for an on-top contact of a bipyridine molecule on a Au (111) surface and as 
parameter for the variation of the surface geometry we used the coordination number Nc (only 



including other Au atoms) for the surface atom directly connected to the molecule (see Fig. CJi 
for three examples). It was also found that the variation of the conductance of such bipyridine 
junctions with Nc is entirely controlled by the position of the transmission peak due to the LUMO 
(Fig. nj)) and that the effect of the interference with other orbitals as well as through- vacuum 
tunneling are negligible. 

But what is governing the LUMO position? We investigate the variation of the energetic 
position of the lowest-lying molecular orbital (denoted MOl in the following text) with respect 
to the metal's Fermi level, in dependence on the distance between the surface and the molecule 
d (Fig. nt). Since MOl is ^ 10 eV below the lowest-lying Au valence states, its position must 
be exclusively guided by rigid potential shifts. Due to the quantitative correspondence of the 
sequence of the LUMO peaks in the transmission functions in Fig. QJ) to the energies of MOl 
at do (Fig. nt), hybridisation effects must be independent of Nc- For large d, MOl rests at 
an energetic position which exactly corresponds to the one it would hold due to vacuum level 
alignment (where the surface structures differ from each other depending on their respective 
work functions), which means that both subsystems are entirely charge neutral. 

Now we want to address the question how the rigid potential shifts that MOl experiences 
at the bonding distance do are related to charge being transferred between the molecule and 
the surface. For this purpose, in DFT the concept of fractional charges can be introduced. ^H] 
This makes it possible to determine the ground state electron density and electronic eigenvalue 
spectrum for a (albeit only finite) system with fractions of electrons removed or added when 
compared to the total charge of all the nucleis. FigOJi shows MOl for isolated bipyridine cations 
in the region where 0.3 to 0.46 electrons have been removed from the molecule. The different 
lines for different values of Nc result because the molecules vacuum level has to be aligned with 
the metals workfunction, which differs for the three cases. If Fig. QJi is now compared with 
the MOl positions at do? partial charges on the molecules can be attributed to different surface 
structures. 

3.2. Biphenyl dithiolate 

For dithiolates there is an on going discussion whether the sulfur atoms of the molecule are 
attached to the Au substrate in a hollow (H) or on top (T) configuration (Fig. Ett). Although 
it is well known, that the hollow configuration is thermodynamically more stable, there are also 
arguments about kinetic effects during the formation of the junction, which would favour an on 
top geometry. 

In a comparison of the transmission functions (Fig. Eb) we find that although the double 
peak structure directly below the Fermi energy is qualitatively similar in both cases, its energetic 
distance from Ejr varies by ^ 0.7 eV. As a direct consequence the conductance for T where 
lies at the shoulder of the peak, is 0.44 Go, whereas for H crossing just the tails it goes down to 
0.03 Go. Inspecting the difference in MOl shifts in Fig. Efc we find that at the bonding distance 
do, it is only ^ 0.3 eV, which indicates that the influence of hybridisation on the position of the 
peaks in the transmission function is quantitatively different for the two geometries. We also 
note that BPDT differs from bipyridine in that charge transfer occurs in the opposite direction. 
The amount of charge being moved at the bonding distance is by almost a factor of two smaller 
for BPDT, which is reflected in much smaller MOl energy shifts in dependence on d. 

The most interesting feature in our study of BPDT is the difference in conductance between 
H and T configurations, which is caused by the energetic variation of the positions of the 
transmission functions in Fig. [^b and the related different lengths where the maxima occur in 
Fig. with respect to do. When d is reduced, the amount of charge transfer increases, which 
makes MOl go up with respect to Ep. At the extremal point the molecule is so close to the 
surface that fractions of electrons are pushed out of the contact region, which is equivalent to a 
net loss of charge. For both, H and T, the position of the maximum of the MOl vs. d functions 
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Figure 2. (Color online) a) Contact geometry for a BPDT molecule connected to a triangle 
of Au atoms on a Au (111) surface. The small spheres denote the positions of the sulfur anchor 
atoms, where H stands for hollow and T for an on top configuration, b) Transmission functions 
and c) energy shifts of MOl for the two configurations in a), where the distance axis is different 
for T and H in absolute values, since it is referring to the respective do in both cases. The inset 
of c) gives the binding energies in dependence on the distance d between the molecule and the 
surface. All energies are given relative to the Fermi energy of the Au surfaces. 

in Fig.Et is at d ^ 2.4. The two bonding configurations only differ in the sense that this distance 
also equals do for T, whereas it is already in the repulsive regime for H (see for comparison also 
the distance dependent binding energies in the inset of Fig. Efc). Therefore it can be concluded 
that the large differences in conductance for the two geometries is an indirect consequence of 
their different bonding behaviours. 

4. Electronegativity 

It is useful to look at our results from the perspective of a quantum chemical approach for 
calculating charge transfer, which is often used in the context of chemical reactions in organic or 
biochemistry The key quantity in this scheme is the electronegativity /x, and supplemented 
with the so called hardness z/, it can be used for computing the total charge transfer between 
any two subsystems A and B, when they are brought into contact, as 

A7V= (/iB-MA)/2(^A + ^B). (1) 

This concept has the appeal that no detailed information about single electron eigenenergies or 
hybridisation of orbitals needs to be analysed or even computed for the coupled system, and 
nevertheless quantitative estimates for the charge transfer can be made. Additionally it includes 
the relaxation of all electrons in the subsystems due to electron gain or loss as a consequence of 
calculating the ionisation potential and electron affinity from total energies of charged molecular 
fragments. As has been pointed out early ^H], however, the screening effects that act on the 
electrons on system B due to the presence of A and vice versa have been neglected in its 
derivation. As a consequence Eq. Q] does not depend on the distance d between the molecule 



and the surface, which is counterintuitive and in contradiction to our results. Recently, a 
semiempirical scheme has been used ^Zj for a qualitative discussion of Fermi level pinning at 
metal-semiconductor interfaces and its dependence on the atomic structure of the interface. 
There a modified version of Eq. d which was corrected for its shape and distance dependency, 
has been derived, 

AN= MB^iM .2) 

2(z^A + J^b) - "^Jab + uaJaa + ubJbb 

In this expression Jab was a parameter for the electrostatic attraction between atoms from 
different sides of the interface, Jaa and Jbb denoted the repulsion between atoms within the 
same plane parallel to the interface and all three quantities were explicitly distance dependent. 
The number of next nearest neighbours within a plane was denoted ua and ub for the two 
subsystems, respectively. If we now compare Eq. [21 with our results for the surface structure 
dependence of the charge transfer for the bipyridine system in Fig. we find that there is 
qualitative agreement in the sense that a small value of our coordination number Nc corresponds 
to a large ua in Eq. El and both leads to an increased amount of transferred charge AA^. 

5. Summary 

In summary, we discussed the level alignment in molecular nanoj unctions for different Au 
surface structures in the case of bipyridine, and for different bonding configurations for BPDT. 
We also compared our findings with the more traditional quantum chemical method of using 
electronegativities for characterising charge transfer, which in its simplest form is based on 
DFT total energy differences and not on single particle electronic eigenvalues. On a qualitative 
level the framework agrees well with our results for bipyridine in dependence on the Au surface 
structure. 
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